Sequential parameter estimation for stochastic systems

The quality of the prediction of dynamical system evolution is determined by the accuracy to which initial conditions and forcing are known. Availability of future observations permits reducing the effects of errors in assessment the external model parameters by means of a filtering algorithm. Usually, uncertainties in specifying internal model parameters describing the inner system dynamics are neglected. Since they are characterized by strongly nonGaussian distributions (parameters are positive, as a rule), traditional Kalman filtering schemes are badly suited to reducing the contribution of this type of uncertainties to the forecast errors. An extension of the Sequential Importance Resampling filter (SIR) is proposed to this aim. The filter is verified against the Ensemble Kalman filter (EnKF) in application to the stochastic Lorenz system. It is shown that the SIR is capable of estimating the system parameters and to predict the evolution of the system with a remarkably better accuracy than the EnKF. This highlights a severe drawback of any Kalman filtering scheme: due to utilizing only first two statistical moments in the analysis step it is unable to deal with probability density functions badly approximated by the normal distribution.